The Anjials of Applied Probability 
2011, Vol. 21, No. 3, 908-950 
DOI: 10.1214/10-AAP711 

(c) Institute of Mathematical Statistics, 2011 



(N 



O: 



(N 



>< 



EFFICIENT SIMULATION OF NONLINEAR PARABOLIC SPDES 
WITH ADDITIVE NOISE 



O " By Arnulf Jentzen^'^, Peter Kloeden^ and Georg Winkel 

Princeton University, Johann Wolfgang Goethe University 
O ■ and Johann Wolfgang Goethe University 



Recently, in a paper by Jentzen and Kloeden [Proc. R. Soc. Lond. 
Ser. A Math. Phys. Eng. Sci. 465 (2009) 649-667], a new method 
for simulating nearly linear stochastic partial differential equations 
(SPDEs) with additive noise has been introduced. The key idea was 
to use suitable linear functionals of the noise process in the numerical 
scheme which allow a higher approximation order to be obtained. 
Following this approach, a new simplified version of the scheme in 

■ the above named reference is proposed and analyzed in this article, 
"j^ , The main advantage of the convergence result given here is the higher 

convergence order for nonlinear parabolic SPDEs with additive noise, 
although the used numerical scheme is very simple to simulate and 
implement. 

■ 1. Introduction. In this article, the numerical approximation of nonlin- 
^ I ear parabolic stochastic partial differential equations (SPDEs) is considered. 

■ Following the idea in [10] for somewhat linear SPDEs, a new numerical 
. method for simulating nonlinear SPDEs with additive noise is proposed and 
\ analyzed in this article. The main advantage of the convergence result in 

this article is the higher convergence order for nonlinear parabolic SPDEs 
with additive noise in comparison to convergence results of classical schemes 
such as the linear implicit Euler scheme. Nevertheless, the here presented 
scheme is very simple to simulate and implement. 

More precisely, let T G (0, oo) be a real number, let (0, j^, P) be a prob- 
ability space and let iJ = L^((0, 1),M) be the M-Hilbert space of equiv- 
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alence classes of square integrable functions from (0, 1) to M. Moreover, 
let / : [0, 1] X M — )• M be a smooth function with bounded partial deriva- 
tives, let ^ : [0, 1] — )• M with ,^(0) = ^(1) = be a smooth function and let 
■.[0,T]xn^ H bea standard Q-Wiener process with a trace class oper- 
ator Q-.H^H (see, e.g.. Definition 2.1.9 in [13]). It is a classical result (see, 
e.g., Proposition 2.1.5 in [13]) that the covariance operator Q-.H^Hol the 
Wiener process : [0, T] x O — )• if has an orthonormal basis gj ^ H, j ^ N, 
of eigenfunctions with summable eigenvalues fij £ [0,oo), j G N. In order to 
have a more concrete example, we consider the choice gj{x) = \/2sin(j7r2;) 
and = cj~^'^~^^^ for all x G (0, 1) and all j G N with some c G [0,oo) and 
some arbitrarily small r G (0, oo) in the following and refer to Section 2 for 
our general setting. Then we consider the SPDE 



,2 



^Xt + f{x,Xt) 



dt + dw9 



t 1 



dXt 

(1) 

Xt(0)=Xt(l) = 0, Xo = e, 

for X G (0, 1) and t G [0,T]. Under the assumptions above, the SPDE (1) has 
a unique mild solution. Specifically, there exists an up to indistinguishability 
unique stochastic process X : [0, T] x — )• with continuous sample paths 
which satisfies 

(2) Xt = e^'i+ f e^^'-'^F{Xs)ds+ f e'^^'-'Uw^ , P-a.s. 

JO JO 

for all i G [0,r] where A:D{A) d H ^ H is the Laplacian with Dirichlet 
boundary conditions on (0, 1) and where F:H^H is the Nemytskii oper- 
ator {F{v)){x) := f{x,v{x)) for all x G (0, 1) and all v £ H. 

Then our goal is to solve the strong approximation problem of the SPDE (1) 
More precisely, we want to compute a J-"/;S(i?)-measurable numerical ap- 
proximation Y -.i}^ H such that 

1/2 

(3) (E 



(k \XT{x)-Y{x)\^dx ^ 



< e 



holds for a given precision e > with the least possible computational ef- 
fort (number of computational operations and independent standard normal 
random variables needed to compute y :f] — )• H). A computational opera- 
tion is here an arithmetical operation (addition, subtraction, multiplica- 
tion, division), a trigonometrical operation (sine, cosine) or an evaluation of 
/ : (0, 1) x M ^ M or the exponential function. 

In order to be able to calculate such a numerical approximation on a com- 
puter, both the time interval [0,T] and the infinite-dimensional M-Hilbert 
space H = L^((0, 1), M) have to be discretized. While for temporal discretiza- 
tions the linear implicit Euler scheme is often used, spatial discretizations are 
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usually achieved with finite elements, finite differences and spectral Galerkin 
methods. For instance, the linear implicit Euler scheme combined with spec- 
tral Galerkin methods which we denote by J-"/S(-ff)-measurable mappings 
Z^:n^H,nG{0,l,...,N^},NGN:={l,2,.. .}, is given by := Pn{0 
and 

^5., := (l-^A 



-'n+l 

(4) 



for every n S {0, 1, . . . , N'^ — 1} and every N £N where the bounded linear 
operators P/v :H ^ H, G N, are given by 

TV 1 
(5) (P/v(f))(x) := y ^2sin(re7rx) / sm{n7r s)v (s) ds 

n=l 

for all x S (0, 1), V £ H and all G N. Note that the infinite-dimensional 
M-Hilbert space H is projected down to the A^-dimensional M-Hilbert space 
Pi\f{H) for the spatial discretization and the time interval [0,T] is divided 
into N"^ subintervals, that is, N'^ time steps are used, for the temporal dis- 
cretization in the scheme Z^ , n G {0, 1, . . . , A^^}, above for iV G N. The exact 
solution X : [0,r] X — )• of the SPDE (1) enjoys at least twice the reg- 
ularity in space than in time and therefore, the quadratic number of time 
steps is used in the scheme (4) above (see also Walsh [15] for details). 

We now review how efficiently the numerical method (4) solves the strong 
approximation problem (3) of the SPDE (1). Standard results in the litera- 
ture (see, e.g.. Theorem 2.1 in Hausenblas [7]) yield the existence of a real 
number C > such that 



(6) E 



\Xt{x) - Z^^{x)\' dx 







1/2 



holds for all G N. Since Pj\f(H) is A^-dimensional and since A'^^ time steps 
are used in (4), 0(A^'^ log(A^)) computational operations and independent 
standard normal random variables are needed to compute Z^2 for A^ G N. 
The log term in 0{N'^ iog{N)) for A^ G N arises due to computing the nonlin- 
earity with fast Fourier transform (aliasing errors are neglected here). Com- 
bining the computational effort 0(A^'^ log(A^)) and the estimate (6) shows 
that the linear implicit Euler scheme needs about 0{£~^) computational 
operations and independent standard normal random variables to achieve 
a precision of size e > in the sense of (3). In fact, we have demonstrated 
that the linear implicit Euler scheme method (4) needs 0{e^^^^^^) compu- 
tational operations and random variables to solve (3) for every arbitrarily 
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small 6 G (0,00) but for simplicity we write about 0(e~^) computational 
operations and random variables here and below. 

Recently, in [10] , a new numerical method for simulating somewhat linear 
SPDEs with additive noise has been introduced. The key idea in [10] is to 
use suitable linear functionals of the noise process in the numerical scheme 
which allows a higher approximation order to be obtained. In this paper, 
we extend this idea to the case of nonlinear SPDEs of the form (1). More 
precisely, we introduce the following numerical scheme which is a simplified 
version of the scheme considered in [10]. Let : $7 — )• if, n S {0, 1, . . . , N}, 

G N, be J^/i3(-ff)-measurable mappings given by Y^^ := Pn{0 ^^'^ 



P-a.s. for every n G {0, 1, . . . , — 1} and every N £N. Note that the infinite- 
dimensional M-Hilbert space H is projected down to the A^-dimensional M- 
Hilbert space Pm{H) for the spatial discretization and the time interval 
[0,r] is divided into N subintervals, that is, N time steps are used, for 
the temporal discretization in the scheme Y^ , n G {0, 1, . . . , A^}, above for 
iVG N. 

We now illustrate the main result of this article (Theorem 1) and show 
how efficiently the method (7) solves the strong approximation problem (3) 
of the SPDE (1). Theorem 1 shows the existence of real numbers > 0, 
6 G (0, 1), such that 



holds for all G N and all arbitrarily small b G (0, 1). The stochastic inte- 
grals 



for n G {0, 1, . . . , A^} and A^ G N in (7) provide more information about the 
exact solution and this allows us to obtain the estimate (8) although only 
A^ time steps (instead of A'^^ time steps in the case of the linear implicit Eu- 
ler scheme) are used in (7). Nevertheless, since the stochastic integrals (9) 
in (7) depend linearly on the Wiener process ^ [0, T] x — )• ii, they 
are again normally distributed and hence easy to simulate. More precisely, 
since P/v(i?) is A^-dimensional and since A^ time steps are used in (7), 
0(A^^ log(A^)) computational operations and independent standard normal 
random variables are needed to compute Y^ for A^ G N. The log term in 



(7) 






(9) 
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0{N^log{N)) for N also arises due to computing the nonlinearity with 
fast Fourier transform (ahasing errors are neglected here). Combining the 
computational effort 0(A^^ log(A^)) and the estimate (8) shows that the nu- 
merical scheme (7) needs about 0(e~^) computational operations and in- 
dependent standard normal random variables to achieve a precision of size 
e > in the sense of (3). 

The estimates (6) and (8) are both asymptotic results since there is no 
information about the size of the corresponding error constants. In particu- 
lar, the error constants Cs £ (0,oo), 6 G (0, 1), in (8) could be much bigger 
than in (6). Therefore, from a practical point of view, one may ask whether 
the numerical method (7) solves the strong approximation problem (3) more 
efficiently than the linear implicit Euler scheme (4) for a given example of 
the form (1) and a given concrete e > 0. In order to analyze this question, 
we compare both methods in the case of a simple reaction diffusion SPDE 
of the form (1) (see Section 4.1 for details) and assume that the strong 
approximation problem (3) should be solved with the precision e = g^g. In 
that example, it turns out that the linear implicit Euler scheme precisely 
needs 2^^ = 2,097,152 independent standard normal random variables while 
the numerical method (7) precisely needs 2^^ = 65,536 independent standard 
normal random variables to achieve an approximation error of size e = ^ 
(see Tables 1 and 2 in Section 4.1). We also emphasize that the numerical 
scheme (7) is very simple to implement and refer to Figure 2 for a short 
MATLAB code. 

Having illustrated the main result of this article, we now sketch the key 
idea in the proof of Theorem 1. The main difficulty was to estimate the 
discretization error for nonlinear F. In that case, the main problem was to 
establish estimates of the form 

^-1 -.(n+l)T/Af 



(10) e^'^^-'\F{X,)-F{X^T/N))ds <Cs-N^'~'^ 

for all G N and all 5 G (0, 1) where Cs G (0,oo), 6 G (0, 1), are appropriate 

constants and where we write ||^||i,2(f7;_f/) '■= (^[Jo 1^(^)1"^ dx])^^"^ G [0,oo] for 
every J^/;S(-ff)-measurable mapping Y li}^ H for simplicity. The smooth- 
ness of the Nemytskii operator F on an appropriate subspace V C H shows 
that it remains to estimate 



(11) 



N-1 

E 

n=0 



{n+l)T/N 



nT/N 



'^F'{X^t/n){Xs-X, 



nT/N 



ds 



L^{Q;H) 



< Cs • N^^~^^ 



for all G N and all 6 G (0, 1). In [10], the linear operators F'{v) for v G 
H and A:D{A) C H ^ H are assumed to commute in some sense which 
is fulfilled in the case of linear F such as F(v) = v, v £ H, but excludes 
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nonlinear Nemytskii operators such as F{v) = (jq^gry, v ^ H (see Assumption 
2.4 in [10] for details). Under this commutativity condition, (11) can easily 
be established by using the smoothing effect of the semigroup e^*, t G [0,r] 
(see Section 5.b.i in [10]). Instead of this condition, our key assumption on 
the nonlinearity is an appropriate estimate on the adjoint operators of the 
Frechet derivative operators of F [see (13)]. Since in our examples F is a 
(nonlinear) Nemytskii operator, the derivative operators F'{v), v £V, are 
self-adjoint and hence, it can easily be seen that this assumption is fulfilled 
[see (17) in Section 4 for details]. Moreover, this assumption enables use to 
show (11) and hence (10) [see Section 6.1.1 and particularly estimate (31)]. 
We also mention that the difficulty to estimate (10) can be avoided by using 
a more complicated scheme with a second linear functional (see Section 6.4 



Finally, we would like to point out limitations of the here presented numer- 
ical method. The following assumption is essential to apply our algorithm. 
The eigenfunctions of the dominating linear operator and of the covariance 
operator of the driving additive noise process of the SPDE must coincide 
and must be known explicitly. 

The rest of this article is organized as follows. The basic setting and the 
assumptions that we use (including our key assumption on the adjoint of the 
Frechet derivative of the nonlinearity) are presented in Section 2. The new 
numerical scheme and its convergence theorem which is the main result of 
this article are given in Section 3. This result is illustrated with some exam- 
ples and some numerical simulations in Section 4. Although our setting in 
Section 2 uses the standard global Lipschitz assumption on the nonlinearity 
of the SPDE, we demonstrate the efficiency of our method numerically for a 
SPDE with a cubic nonglobally Lipschitz nonlinearity in Section 5. Proofs 
are postponed to the final section. 

2. Setting and assumptions. Fix T G (0,oo) and let (r2,-F, P) be a prob- 
ability space with a normal filtration {J~'t)te[o,T] which means J-t+ = J-'t 
for ah t G [0,r) and {A £ J"|P[A] = 0} C J"o (see, e.g.. Definition 2.1.11 
in [13]). In addition, let {V, \\ ■ ||y) be a separable R-Banach space and let 
{H, {■,-)h, II ■ \\h) be a separable M-Hilbert space with V C H continuously. 
The following assumptions will be used. 

Assumption 1 (Linear operator A). Let {Xn)n&N C. (0, oo) be an in- 
creasing sequence of real numbers and let (e„)„gN C be an orthonormal 
basis of H. Assume that the linear operator A : D{A) C H ^ H is given by 



in [11]). 



oo 




n=l 



for all v£D{A) with D{A) = {w £ H\Y,n=i\^n\^\{('n,w)H\'^ < oo}. 
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Let D{{-AY) with \\v\\d{{^aY) = \\{-AYv\\h for v G D{{-AY) and r G M 
denote the domains of fractional powers of the hnear operator —A (see, e.g., 
Section 3.7 in [14]). 

Assumption 2 (Nonhnearity F). Assume _D((— A)^/^) C V continuously 
and let F -.V ^ V be a twice continuously Frechet differentiable mapping 
with 

||-F'(f )w||/f < c||i(;||/f, 

(12) 

\\F'{v)\\Liy)<C, ||F"(t;)||^(„(^)<C, 

(13) ll(^'(^i))*llL(D((-A)i/2)) < C(l + lk||^((_^)i/2)) 

for every v,w €zV and every u G D((— A)^/^) where c G [0, oo) is a given real 
number. 

By definition F'(v) G L{V) is a bounded linear mapping from V to V for 
every v gV. Due to the first condition in (12), we also have that F'{v) G 
L[H) is a bounded linear mapping from H to H for every t; G In that 
sense, the adjoint operator {F'{v))* G L{H) given by 

((F'(^;))*n,i(;)H = (n,F'(z;)w;)H 

for all u,w £ H is well defined for every v €V. Due to (13), the opera- 
tor {F'{v))* G L{H) is also a bounded linear mapping from D{{-AY/^) to 
D{{-AY/^) for every u G D{{-Ay/^). 

Assumption 3 (Stochastic process O). Let O: [0,T] x ^ D((-^)^) 
be a centered and adapted stochastic process with continuous sample paths 
such that Ot2 - e^(*2-ti)(9j^ independent of Tt^ for aU < ti < ^2 < T and 
such that 



E 



'-0<t<T 



sup \\{-AyOt\\jj + sup iit2-hy^'E[\\Ot,-OtXv])<oo 



0<ti<t2<T 



holds where 7 G [^,1) and 9 G (0, ^] are given real numbers. 

Assumption 4 (Initial value ^). Let D{A) be a To/B{D{A))- 

measurable mapping with E[||A^|||^] < 00. 

These assumptions suffice to ensure the existence of a unique solution of 
the SPDE (14). 

Lemma 1 (Existence of the solution). Let Assumptions I-4 be fulfilled. 
Then there exists a unique adapted stochastic process X : [0, T] x $7 — D{{—A)"') 
with continuous sample paths which fulfills 

(14) Xt{uj) = e^'^{u)+ [\^('-'^F{Xs{uj))ds + Ot{u:) 

Jo 
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for all t G [0, T] and alluj^n. Moreover, X:[0,T]xn^ D{{-A)'^) satisfies 
E[supo<t<Tll(-^r^t|||^]<oo. 

The proof of Lemma 1 is given in Section 6. Some examples satisfying 
Assumptions 1-4 are presented in Section 4. 

3. Numerical scheme and main result. For numerical approximations of 
the SPDE (14), we have to discretize both the time interval [0,r] and the 
M-Hilbert space H. To this end, we use projections P/v :H ^ H given by 
■> 'v) H s-n for every v E H, N (z ^ and finite-dimensional 
M-Hilbert spaces Hn C H given by ■= Pn{H) for every iV G N. Fi- 
nally, we define J-"/i3(ifAr)-measurable mappings 

{0, 1, . . . , M} and N, M G N by Y^'^\lo) := Pn{^{uj)) + Pjv(Oo(w)) and by 



VN,M 



(15) 



+ -fW(0(m+i)T/A/('^) - e^^/^^0„T/A/(^)) 



for every m G {0, 1, . . . , M — 1}, N,M G N and every a; G In many exam- 
ples, this scheme is as easy to simulate as the classical linear implicit Euler 
scheme. We refer to Section 4 for a detailed description of the implementa- 
tion of our numerical scheme including a short MATLAB code. 

Theorem 1. Let Assumptions 1~4 be fulfilled. Then there is a real num- 
ber C > such that 

(16) (E[iix.„,,„ - Y,;-'-iay^ Hx^^ ^i^^) 

holds for every m G {0, 1, . . . , M} and every N,M e where {Xn)n£N C 
(0,oo) is given in Assumption 1 where 7 G [^,1) and 6 G (0,^] are given 
in Assumption 3 where X:[0,T] x ^ D{[—A)'^) is the solution of the 
SPDE (14) and where Y^'^'^ : n Hn, m G {0, 1, . . .,M}, N,M G N, is 
given by (15). 

Here and below log is the natural logarithm. While the expression jjj^ 
for G N in (16) arises due to discretizing the infinite-dimensional R-Hilbert 

space H, the expression ^^^j^^e^^^^ for M G N arises due to discretizing the 
time interval [0,r]. We would like to remark that the logarithmic term 
in for M G N can be avoided by assuming F{D{{-Ay/^)) C 

D{{—AY) and an appropriate linear growth condition on F for some e > 0. 
Although this condition is fulfilled in our examples below, we use this loga- 
rithmic term in Theorem 1 here in order to formulate Assumption 2 in our 
abstract setting as simple as possible. 
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A similar result could be obtained for SPDEs of the form (14) but with a 
time dependent nonlinearity F. However, we omit the time dependency of 
the nonlinearity here for simplicity. 

4. Examples. Let i7 = L2((0, 1),R) be the M-Hilbert space of equivalence 
classes of ;B((0, l))/B(M)-measurable and square integrable functions from 
(0,1) to M with the scalar product and the norm given by 

1/2 



{v,w)h 



v{s)w{s) ds, 



\v\\h 



\v{s)\Us 



for every v,w £ H. In addition, let V = C([0, 1],M) be the M-Banach space 
of continuous functions from [0,1] to M equipped with the norm \\v\\v = 
supo<a;<i|f (x)| for every v &V. 

Let K £ (0,oo) be a given positive real number and let {Xn)neN C (0,oo) 
and {en)neN C H he given by 



A. 



2 2 

KTi vr , 



e„(x) := \/2sin(n7rx) 



for every x S (0, 1) and every n G N. Hence, the linear operator A : D{A) C 
H ^ H reduces to the Laplacian with Dirichlet boundary conditions on the 
interval (0,1) times the constant k £ (0,oo) (see, e.g.. Section 3.8.1 in Sell 
and You [14]). In particular, £>((— A)^/^) reduces to the M-Sobolev space 
i/(J((0, 1),M) equipped with the norm 

;-^)V2, 



lD({-A)i/2) 



U H 



\n=l J 



1/2 



1 \ 1/2 

\'\i(x)^ dx 



for all u G L»((-A)V2). (See Sell and You [14] for more information about 
this space.) 

Furthermore, let / : [0, 1] x M — t- M be a twice continuously differentiable 
function with the bounded partial derivatives 

{x,y) 



dl 
dy 



dx dy 



<K 



for all x € [0, 1] and all y G R with an arbitrary constant K G [0, oo). Then 
the Nemytskii operator F :¥ given by (F(v))(x) = f{x, v{x)) for every 
X G [0, 1] and every v gV satisfies Assumption 2. To see this note that 

^^(x,n(x)) -vix), 



F\u){v) 
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) & V . Therefore, we have 



holds for all u,v^w 

= \\{-Ar'F'{u)v\\l 

' d_ 

dx 



-- K 



= K 



< 2k 



(x, u{x)) + (^^^ (x, u{x)) ■ v'{x] 



^ ^ -^-^{x,u{x)) '^v{x) dx 



dx 



and 



+ 2k 



d f \ ^ 
— ] (x,u{x)) ■ v' (x) dx 
dyj 



<2k\\v\\1 



d 



dx 



^^{x,u{x)) dx 



+ 2kK^ [ \v'{x)\^dx 
Jo 

<4|bllVA)V.)/ I (^) 



dx 



+ 4|I^IIL((-A)V.) I 

+ 2i^^lbllL((_A)i/2) 

Hence, we obtain 



\dxdy J 



dx 



for all u,veD{{-A)^/^) 

+ 2i^^|b||^((_^)l/2) 



- 6-?^^||'v||^((_^)i/2) +4K^K ^lbll^((_A)i/2)ll'"llD((-yl)i/2) 
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and 



||(F'(n))*t'||£,((_^)i/2) 



|2 +47<'2k-1||7;||2 ||,,||2 



(17) 



< y/QK\\ 



< K\\v\ 



^llD((--yl)i/2) 



lz)((-A)l/2)||U'|lz)((_yl)l/2) 



lD{(-A)i/2)(3 + ^'"ll^^llD((-yl)l/2)J 

for all u,v £ D{{—A)^^'^). This shows that F indeed satisfies Assumption 2 
with c = 3 A'. 

Let (6n)neN C M be a sequence of real numbers with Yl'^=i n'^ \bn\'^ < oo iov 
some arbitrarily small e G (0,oo). Lemma 4.3 in [1] then gives the existence 
of an up to indistinguishability unique stochastic process O: [0,T] x — )• ^ 
which satisfies Assumption 3 for = ^ and 7 = ^ and which satisfies 



-1/2 I 



Ot 



n=l 



1 



for all t £ [0,T] where the /3" : [0,r] x R, n£ N, are independent stan- 
dard Brownian motions with respect to a given normal filtration {J't)te[o,T]- 
Moreover, the mapping ^ : $7 — )• given by 

(C(a;))(x) = -ijsin(7rx) + ^sin(37rx) 

for all w G and all x G (0, 1) obviously satisfies Assumption 4. 
In view of the above choice, the SPDE (14) reduces to 

;Xt + f{x,Xt) dt + BdWt. 



(18) 



Xt{0) = Xt{l)=0, 

, , sinfvrx) 3\/2 
Xo(x) = — ^ + sin(37rx) 



' 5 

for X G [0, 1] and t G [0, T\ where the linear operator B iH ^ H is given by 

oo 

Bv = '^bn{en,v)en 

n=l 

for ell V G H and where {Wt)te[o,T] is a cylindrical /-Wiener process on H. 

Since Assumption 3 is fulfilled for = ^ and 7 = ^ , Theorem 1 shows the 
existence of a real number C > 0, such that 



(19) E 



\X^nT/Mi^)-Y.^''''ix)\'dx 



1/2 



< 



fl (l + log(M)) \ 
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holds for all n G {0, 1, . . . , M} and all N, M G N. While the expression for 
A'^ G N in (19) corresponds to the spatial discretization error, the expression 

^'^"'"'m^^'*^ M G N in (19) corresponds to the temporal discretization error. 
Since these error terms are nearly of the same size, we choose M = N and 
consider the numerical approximations Y^'^ : — t- Hj\[, n G {0, 1, . . . , A^}, 
G N, in the following. Due to (19), we obtain the existence of real numbers 
> 0, (5 G (0, 1), such that 



(20) E 







\xT{x)-Y^'''{x)\^dx 



1/2 



holds for all G N and all arbitrarily small 6 G (0, 1). 

In order to describe the implementation of the numerical scheme (15) in 
this example, we use the J^/;B(]R)-measurable 

by Y^^fiu) := {en, Y^'^''{oj))h for ah n G {1, 2, . . . , A^}, m G {0, 1, . . . , M} 
and all A^, Af G N. The numerical scheme (15) for the SPDE (18) with M 
N tl 
and 



A^ then reduces to Y-^q^ = ^, Y^^^^ = 0, ^3^'^ = f , ^40'^ = ^5^'^ = • • • = 



Oi • tt\/2k 



^2,n+l-e lJ^2,n + ]^ (^2 , ( ))Hj 



(21) + ^ r- 



XN,n 



hN ■ A^7rV2K 

for all n G {0, 1, . . . , A^ - 1} and all A^ G N where the J'/B(M)-measurable 
mappings Xn,m : ^ M for n G {1, 2, . . . , A^}, m G {0, 1, . . . , A^- 1} and A^ G N 
are independent standard normal random variables. Since 0(A^^ log(A^)) 
computational operations and independent standard normal random vari- 
ables (computational effort) are needed to compute the numerical solu- 
tion Yj^ ' given by (21) for A^ G N, it fohows that Y^''' converges with 
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order ^ with respect to the computational effort to the exact solution 
X:[0,T]xn^ L>((-A)i/2) of the SPDE (18) in the sense of (20). We re- 
mark that the log term in the computational effort 0{N'^ log(A^)) for G N 
arises if one computes the nonlinearity in (21) with fast Fourier transform 
(see Figure 2 for details). 

In order to compare the new numerical scheme (21) with classical schemes, 
we consider the well-known linear implicit Euler scheme combined with 
spectral Galerkin methods applied to the SPDE (18). The linear implicit 
Euler scheme is denoted by J^/i3(i/Ar)-measurable mappings : Q. — )• i^Ar, 
n G {0, 1, . . . , N^}, e N, given by {u:) := P7v(^(a;)) + Piv(Oo(^)) and 



Z^_,,{c.) ■.= {!- ^a) {z^^iu) + I^iP^DiZ^iu)) 



(22) 



Ar2 



for every n G {0, 1, . . . , A^^ — 1} and every N G'N. It has been shown in the 
literature (see, e.g., Walsh [15], Gyongy [4] and Hausenblas [7]) that the 
linear implicit Euler scheme (22) and other classical numerical schemes such 
as the linear implicit Crank-Nicolson scheme combined with finite elements, 
finite differences and spectral Galerkin methods converge with order | with 
respect to the computational effort. 

The following two numerical examples illustrate the convergence order 
^ of the numerical scheme (21) and the convergence order | of the linear 
implicit Euler scheme (22). 



4.1. A stochastic reaction diffusion equation. In this example, we set 



K- — -!- T— 1 h — " 

~ 100' — J-i "rt — 3.5 



for all n G N and consider / : [0, 1] x 



given 



by f{x,y) = 5(^^ for all x G [0, 1], y G M. The SPDE (18) then reduces to 



dXt 



(23) 



100 dx"^"" ' "(1 + X| 
Xt{0) = Xt{l) = 0, 

Xq{x) = -j= 1 — sin(37rx) 

V 2 5 



dt^BdWu 



for X G [0, 1] and t G [0, 1]. In Figure 1 (see also Tables 1 and 2), we plot the 
root mean square discretization error 



(24) 



E 



\Xt{x)-Y^^''{x)\Ux 



uo 



1/2 
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Computational effort * log{M) Computational effort * log(N) 



Fig. 1. Root mean square approximation error (24) of the numerical scheme (21) and 
root mean square approximation error (25) of the linear implicit Euler scheme (22) applied 
to SPDE (23) versus up to a constant the computational effort. 



Table 1 

Root mean square approximation error (24) ofY^'^ given by (21) 
applied to the SPDE (23) for N G {2^ 2^ . . . , 2^^} 



Numerical 
scheme (21) 


Independent standard 
normal random 
variables N'^ 


Computational 
effort log(iV) 
(up to a constant) 


Root mean square 
approximation 
error (24) 


^22 


16 


22 


0.1864 


V2^.2^ 
^23 


64 


133 


0.0914 


y2^2* 
-^2^1 


256 


710 


0.0417 


T/-2=,2S 
-^25 


1024 


3549 


0.0191 


V2''.2'^ 
-■^28 


4096 


17,035 


0.0091 


2'' 


16,384 


79,496 


0.0045 


V2*.2* 
-^28 


65,536 


363,408 


0.0022 


V2^2^ 


262,144 


1,635,339 


0.0011 


V2 .2 
-I 210 


1,048,576 


7,268,174 


0.0005 


V2",2" 
-•^211 


4,194,304 


31,979,969 


0.0003 



of the numerical scheme (21) versus A^^log(A^) (up to a constant the com- 
putational effort) and the root mean square discretization error 

1/2 



(25) 



E 



\XTix)-ZUx)\^dx 



EFFICIENT SIMULATION OF SPDES 



15 



Table 2 

Root mean square approximation error (25) of Z^2 given by (22) applied 
to the SPDE (23) /or iV e {2\ 2^, . . . , 2^} 



Linear implicit 
Euler 

scheme (22) 


Independent standard 
normal random 
variables 


Computational 
effort iV^ log (AT) 
(up to a constant) 


Root mean square 
approximation 
error (25) 




8 


6 


0.3066 




64 


88 


0.1715 




512 


1064 


0.0837 


^28 


4096 


11,356 


0.0353 


'72^ 
^210 


32,768 


113,565 


0.0135 




262,144 


1,090,226 


0.0058 


^2" 


2,097,152 


10,175,444 


0.0027 



of the linear implicit Euler scheme (22) versus A^^log(A^) (up to a constant 
the computational effort) for different G N. The "expectations" are based 
on 40 independent random realizations and the unknown "exact" solution 
is approximated with a very high accuracy there. 

The short MATLAB code in Figure 2 shows that the solution of SPDE (23) 
can be simulated quite easily with the numerical scheme (21). Figure 3 is the 
result of the matlab code in Figure 2. It shows the solution of the stochastic 
reaction diffusion equation (23) at time t = T = 1 for one sample path c<; G 
approximated with the numerical method (21). 

4.2. A stochastic partial differential equation with a spatially dependent f . 
This time let k = T = 1, bn = — for all n € N and consider / : [0, 1] x 
M M given by f {x,y) = (3.8x^ - 2)y for all x G [0, 1], y G M to obtain the 
SPDE 



dXt 



l^X. + (3.8.^-2)X, 



dt + BdWt, 



N = 1000; T = 1; A = - pi-2 * (l:N).-2 / 100; Y = [l/2,0,3/5,zeros(l,N-3)] ; 
S = sqrt( ( exp(2*T/N*A) - 1 ) ./ A / 2 ) / 3.5 .* (1:N).- -0.55; 
for n=l:N 

y = dst(Y) * sqrt(2) ; 

FY = idst( 5*(l-y) ./(1+ y.-2 ) ) / sqrt(2); 
Y = exp( A * T/N ) .* ( Y + T/N * FY ) + S . * randnCl.N); 
end 

plot( (0:N+1)/(N+1) , [0,dst(Y)*sqrt(2) ,0] , 'k' , 'Llnewldth' , 2 ); 

Fig. 2. MATLAB code for the numerical scheme (21) applied to the SPDE (23). 
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1.4 




Fig. 3. Result of the matlab code m Figure 2: Solution of the stochastic reaction diffu- 
sion equation (23) at t — T — 1 for one sample path to £Q approximated with the numerical 
method (21). 



(26) 



xm=xt{i)=Q, 



Xo(x) 



3\/2 



sm(37rx) 



V2 5 

for X G [0, 1] and t G [0,T]. Here too, the numerical approximation (21) con- 
verges to the exact solution with order ^ with respect to up to a constant 
the computational effort (see Figure 4). Finally, in Figure 5 we illustrate 
how the two different / from examples (23) and (26) affect the evolution of 
the respective solution Xt{oj,x), x G [0,1], for t G {0)^)^)^)1} a-nd one 
sample path lo £^1. 



5. A further numerical example. Although our setting in Section 2 uses 
the standard global Lipschitz assumption on the nonlinearity of the SPDE, 
we demonstrate the efficiency of our method numerically for a SPDE with 
a cubic nonglobally Lipschitz nonlinearity in this section. More formally, we 
consider the SPDE 



(27) 



with 



and 



dXt 



1 

lb 



92 



92 



dt + dW^ 



t ' 



X, 



tl9(0,l)2 



Xq{xi,X2) = sin(7rxi) sin(7rj;2) 
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P 






Numerical scheme (21) 
■• O Orderlines 1/3, 172 





Computational effort * log{N) 

Fig. 4. Root mean square approximation error (24) of the numerical scheme (21) applied 
to SPDE (26) versus up to a constant the computational effort. 



for t,xi,X2 G [0, 1] on the M-Hilbert space H = L2((0, 1)^M) of equivalence 
classes of ;B((0, l)^)/;B(IR)-measurable and square integrable functions from 
(0, 1)^ to M here where (Wj*^)jg[o^i] is a cylindrical Q- Wiener process on H 
with the covariance operator Q-.H^H given by 

{Qv){xi,X2) 

4sin(?i7r2;i) sin(?i7r2;2) 

n,m= 



(n + m)^ 

"1 



sin(n7ryi) sin(m7ry2)i;(yi, y2) dyi dy2 

10 JO 

for all xi,X2 G (0, 1) and all v £ H . Of course, (27) is not included in our 
setting in Section 2. Even worse, it has recently been shown in [9] that 
many numerical methods fail to converge to the solution of a stochastic 
differential equation with super linearly growing coefficients in the strong 
root mean square sense. However, convergence in the pathwise sense often 
holds due to Gyongy's result [3]. Therefore, we plot in Figure 6 the pathwise 
difference 

rl rl \ 1/2 

/ / \Xt{u},Xi,X2) -Y^'^{uj,xi,X2)\'^ dxidx2 
Jo Jo 
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Stochastic reaction diffusion equation (23) SPDE (26) 



Solution t = 
T.5f ■ ■ — 




0.2 0.4 0.6 0.8 1 



X 

Solution t = 0.1 




0,2 0.4 0.6 0.8 



X 

Solution t = 0.3 








0.2 0.4 0.6 0,8 1 



X 

Solution t = 0.6 




0,2 0.4 0.6 0.8 1 

X 

Solution t = 1 

1.5 1 ■ ■ ■ 1 








0,2 0,4 0.6 0,8 1 
X 

Fig. 5. Solution Xt{uj,x), x £ [0,1], of th 
and of the SPDE (26) /or t £ {0, jij, ^, ^, : 
with the numerical method (21). 



Solution t = 




0.2 0.4 0.6 0.8 1 



X 

Solution t = 0.1 




0.2 0.4 0,6 0.8 1 
X 

Solution t = 0.3 




0.2 0.4 0.6 0.8 1 



X 

Solution t = 0.6 




0,2 0,4 0,6 0,8 1 

X 

Solution t = 1 

1.5 1 ■ ■ ■ ■ 1 




0.2 0,4 0.6 0.8 1 
X 



stochastic reaction diffusion equation (23) 
and one sample path a; G SI approximated 



of the exact solution and of the numerical approximation Y^' (w) 

[see (15)] appHed to the SPDE (27) versus up to a constant the computa- 
tional effort N^\og{N) for N £ {2^, 2^, . . . , 2^} and one random cj G It 
turns out that the method (15) converges with order | with respect to the 
computational effort. The linear implicit Euler scheme is known to converge 
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10° 10== 10' lO*' 10" 



Computational effort * log(N) 

Fig. 6. Pathwise approximation error of the numerical scheme (15) applied to SPDE (27) 
versus up to a constant the computational effort for one random a; £ SI. 



in the pathwise sense with order ^ with respect to the computational effort 
to the solution of the SPDE (27). Further pathwise approximation results 
for the SPDE (27) and other SPDEs with nonglobally Lipschitz coefficients 
can be found in [4-6] and [12], for instance. Finally, we plot the solution of 
SPDE (27) for t G {0, ^} and one random w G in Figure 7. 

6. Proofs. The notation 

II^IIl.(w:=(E[||^||^])^/^ G[0,oo] 

is used throughout this section for an M-Banach space iW, \\ ■ \\w)-, a /BiW)- 
measurable mapping Z :0, and a real number p G [1, oo). 

6.1. Proof of Theorem 1. The J-'/fi(//)-measurable mappings : 17 — ?• 
for m G {0, 1, . . . , M} and M G N given by 

(rn-l \ 
^e^('"'^-'='^)F(Xfc,(a;)) 
^ ^ k=o ) 

(28) 

for every m G {0, 1, . . . , M}, G and M G N are used throughout this 
proof. Here and below h is the time stepsize h = Hm = with M G N. This 
proof is divided into three parts. In the first part (see Section 6.1.1), we 
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Solution t = 




Fig. 7. Solution Xt{Lo, xi,X2), xi,X2 G [0,1], of the SPDE (27) for t £ {0, ^} and one 
random oj £fl approximated with the numerical method (15). 



estimate 

II V "V^A^^'ll 

W^mh — im llL2(n;_ff) 

for every m E {0, 1, . . . , M} and every M E N which corresponds to the tem- 
poral discretization error. In the second part (see Section 6.1.2), we estimate 

\\^m - PN{ym)\\L-^{Q.;H) 
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for every m G {0, 1, . . . , M} and every A^, M G N which corresponds to the 
spatial discretization error. Finally, we estimate 



m ) 



N.M\ 



for every m E {0, 1, . . . , M} and every N,M in the third part (see Sec- 
tion 6.1.3). Combining these three parts will then yield the desired assertion 
via Gronwall's lemma as we will see below. 

Before we begin with the first part, we introduce a universal constant 
i? > which is needed throughout this proof. More precisely, let Rg (0,oo) 
be a real number which satisfies 

\\F{Xt)\\L2in;H)<R, 

\\{Xt, -Ot,)- {Xt, - OtML-{n;H) < R\t2 -hi 



1 1 

a7 + (T^ 



+ T + c<R, 



10*2 - OtJ|i4(Q.y) < R\t2 - hf, 
\\^\\L'2(n:D((-A)'y) < ^1 



t|lL4(n;D((-A)i/2)) < R 

for every t,ti,t2 G [0,T] and every v gV where Ai € (0,oo) is given in As- 
sumption 1 where c G [0,oo) is given in Assumption 2 and where 7 G [^,1) 
and 9 G (0, ^] are given in Assumption 3. Indeed, such a real number exists 
due to Assumptions 1-4 and Lemma 4 in Section 6.2. 

6.1.1. Temporal discretization error. Due to (14), we have 

rmh 

^''^'^^-'^F{X,)ds + Om.h 



Xrnh = e^^''i + 





m— 1 



Amh 



A:=0 



{k+l)h 



,A{mh—s) 



kh 



F{Xs)ds + 0. 



mh 



for every m G {0, 1, . . . , A/} and every M G N. From (28), we have 

\\Xmh - Ym\\L'^{Q.;H) 
m-l p{k+l)h 

'kh 



fc=0 



m— 1 



A{mh-s)p^j^^^ (is - /i e^^"'''-''^^F{Xkh) 



.k=0 



L2{n;H) 



22 



A. JENTZEN, P. KLOEDEN AND G. WINKEL 



< 



^-2 j-{k+l)h 



E 



fc=0 



kh 



+ 



mh 



^A{mh—s) 



max(m— l,0)/i 
Ah 



Fix, 



lL2{f7;H) 



k=0 

ds 



+ h\\e -F(^max{m-l,0)h)llL2(n; 



and 



\Xm.h ~ ^ 



Ml 



^mh — J^m \\L^{Q;H) 
"1-2 r(k+l)h 

E 

fc=0 



< 



fc/i 



(m-2 
eMmh-kh)p^X^^^ 
k=0 



+ 



mh 

max(m— l,0)fe 
Ah 



\e^^"''-'^\LiH)\\nXs)\\LHn;H)ds 



+ h\\e ||L{H)l|-^(-^max(m-l,0)h)llL2(n.^) 

for every m E {0, 1, . . . , M} and every M G N. Therefore, we obtain 



\Xmh — Ym \\L'^{n;H) 
"^"2 r.{k+l)h 



< 



E 

A;=0 



^A{mh— s) 



kh 
m-2 



F{Xs)ds 



h[Y,e^^'^''-'"'^F{Xkh) 



.k=0 



+ 2Rh 



L^(n;H) 



< 



"■-2 p(k+l)h 



kh 



E 

k=0 

+ 2Rh 

"1-2 r(k+l)h 

+ E 



A{mh—s) 



iF{Xs)-F{Xkh))ds 



k=0 



^A(mh—s) 



F{Xkh)ds 



kh 



''m-2 



ME 



^A{mh-kh) 



Fix, 



kh 



.k=0 



L2(Q;H) 



and 



-Ml 



l^mh — Ym \\L^{n;H) 



{k+l)h 



< 



m-2 

E 

fc=0 



,A{mh— s) 



iFiXs)-FiXkh + Os-Okh))ds 



kh 



L2{n;H) 
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+ 



"1-2 /•{k+l)h 



E 

k=0 



^A{mh—s) 



(FiXkh + - Okh) - F{Xkh))ds 



m~2 

+ E 

fc=0 



kh 



{k+l)h 



kh 



L'2{n;H) 

+ 2Rh 



L^{n;H) 



for every m G {0, 1, . . . , M} and every M G N. Hence, we obtain 



W^mh — Ym\\L'^{n; 



m-2 
fe=0 



{n;H) 
{k+l)h 



^A{mh—s) I 



+ 



+ 



kh 

^(fc+l)/i 



LiH)\\Fi^s) - F{Xkh + 0s- Okh)\\LHn;H) ds 



E / e^(-'^-^)F'(Xfc,)(0, - Okh) ds 

k=o ■'^^ 

'^-2 ,.{k+l)h pi 



X {Os - Okh, Os - Okh) 
X (1 — r) drds 



L'2{n;H) 



m-2 

+ E 

A;=0 



(fc+l)/i 



gyl(m/i— s) A{mh—kh) 



2 )\#— 1 



)F{Xkh)\\L,^n;H)ds + 2R'M 



and 



I 



Ml 



-'^s - (-'^fch + Os - Okh)\\L2(n;H) ds 

fc=0 

™-2 ,.{k+l)h 



+ 



+ 



1 — n 



fc=0 



m-2 



fe=0 '^'^ 



L^{n;H) 



™-2 />(fc+l)h. 

+ E 



fc=0 



A(mh— s) II 

e llL(/f) 
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X [\\F"{Xkh + r{Os-Okh)) 
Jo 

X {Os - Okh,Os - Okh)\\L2(n;H) drds 



{k+l)h 



m-2 
k=0 '^^'^ 



for every m € {0, 1, . . . , M} and every M G N. Therefore, we have 

"^-2 i-{k+l)h 

<cY, \\[Xs-Os)-[Xkh-Okh)\\L^in;H)ds 

k=0 ■^^^ 



+ 



"1-2 ,.{k+l)h 
k=0 ■^^^ 



L'2{n;H) 



™-2 ^(k+l)h 



k=0 



+ cRY: / / lllio, 



s — Ofc/i|lv'|lL2(Q.K) drds 



k=0 



m-2 



{k+i)h _ ^/^ _ + 



fc=0 
?2 )\/r— 1 



— — r \\F{Xkh)\\LHn;H) ds 

[mh — s) y ' / 



+ 2R^M- 



for every m E {0, 1, . . . ,M} and every M G N due to Lemma 2 below (see 
Section 6.2). Furthermore, we have 



E 



5^ / e^(-''-^)F'(Xfc;,)(0,-e^(^-^'^)0,;,)ds 
fe=0 '^^ 

™-2 r / f{k+l)h 



fc,fc=0 
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and hence 

„(k+l)h 



E 



7 f) /i^/l 

'^-2 r p{k+l)h 
7 — n L -//c/l 



k=0 
m-2 



H. 



m-2 

k,k=0 
k^k 



{k+l)h 



H 



kh 

{k+l)h 
kh 



E 



I H 

for every m G {0, 1, . . . , M} and every M G N. This yields 

™-2 f(k+l)h ^ 
k=0 

™-2 p /.(fc+l)/i 

= E / e^(™'^-^)F'(X,;,)(0, - e^^'-'"'^0,h) ds 

k=0 L 

"^-2 r / (■(k+l)h 
+ 2 V E ( / e^("^'*-")F'(Xfeft) 



2 



k,k=Q 
k<k 



C{k+l)h 

/ e^'^'^^-'^ F' {Xj^ 

Jkh 



H 



and 



E 



"1-2 p{k+l)h 



fc=0 



H. 
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k,k=0 
k<k 



{Os-e''^'-'"'^Okh)ds, 



[■(k+l)h 
Jkh 



for every m € {0, 1, . . . , M} and every M G N. Hence, we obtain 



E 



r.{k+l)h 
'71-2 r f{k+l)h 

fc=o 

"1-2 p , f{k+l)h 

+ 2 V E ( / e^("'''-")F'(Xfc;,) 



fc=0 
m-2 



2 - 

H. 



k,k=0 
k<k 



(E[0.-e^(-^'^)O^JJ-J)d. 



H 



and 



E 



(30) 



™-2 j-(k+l)h 



fc=0 
m-2 



k=0 



{k+l)h 



kh 



HA 
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for every m G {0, 1, . . . , M} and every M S N due to Assumption 3. Com- 
bining (29) and (30) then shows 



\Xmh — Yrn \\L'^{Q,;H) 



{k+l)h 



m.~2 



fc=0 



kh 



\\{Xs - Os) - {Xkh - Okh)\\L^{n;H) ds 



m-2 2 



>fc=0 



kh 



1/2 



"1-2 „(fc+l)/i 



kh 



\L^(n;H) 



ds 



k=0 

+ cR^ \\0s-0kh\\l4^n.v)ds+^ 

k=0 '^^^ fc=0 ''^ 



k=0 



it 7 — ; as 

'kh [mh-s) 



for every m E {0, 1, . . . , M} and every Af G N. Hence, we obtain 

\\Xmh - '^m\\L'^(Q.;H) 



m-2 



k=0 
m-2 

+ E 

A:=0 
/ m 

+ E 



{k+l)h 



kh 



(s — kh) ds 



'(m/i- (A; + l)/i) 



ds 



"1-2 2 

e^(™'^"^)F'(Xfc;,)(0, - e^<^'~'"''>Okh) ds 

kh 



.k=0 



1/2 



m-2 

+ E 

A;=0 



(k+l)h 



^A{mh—s) 



1/2 I 



\LiH) 



X ||(-^)-V2^/(^^,^)((gA(.-fc/^) 
'J^-2 p{k+l)h 
Ikh 



»n-2 

+ cRY, \\Os-Ouh\\l4(a,v)ds + 2R^M-^ 

k=Q 



and 



< ^cRMh'^ 
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m—2 , 
^ 2(m -k-l) 
( "1-2 / f{k+l)h 



1/2 



"1-2 f.{k+l)h 
' ^ 1 1. 1. 



Ikh 

X 



"1-2 Ak+l)h 

+ cRY {R{s-khffds 

k=0 ''^'^ 

for every m S {0, 1, . . . , M} and every M G N. This yields 

\\Xmh — '^m\\L'^{n;H) 
1 1 

\fc=i / 

"*-2 f{k+l)h 

k=o ''^'^ 

(m-2 , r.{k+l)h 
™-2 f(k+l)h 

/ (m/i-(fc + l)/i)^^/2 



+ 
fc=o 



< - 

- 2 



V k=2 J 

+ 2R^M-^ + cR^Mh^^+'^^^ 
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{™-2 ,.{k+l)h ~\ 
k=o ■'^^ ) 

■rn-2 f{k+l)h 

x({e''"-""--f)0/n)|li,,n;«,<i» 
for every m € {0, 1, . . . , M} and every M G N. Hence, we have 

\\Xmh - ym\\L'2{n;H) 

1 1 / f^^ 1 \ 



+ 



. k=0 

m-2 



(m-fc-l)/i 



Ikh 



X 

and 



< ^R'^M-^ + ^R^M-\l + log(M)) + 2R^M-^ + R^M''^^ 

{™^2 r.{k+l)h 
fc=0 •^^'^ 



m-2 ^ 



(m- /c- l)/i 

r-(fc+l)/i 

sup \{w,{-Ar^'^F'{Xuh) 

\\M\h<1 



X 



for every m E {0, 1, . . . , M} and every M G N. This yields 

\\Xmh - Ym\\L'^{Q.;H) 



L2(n;IR) 
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m-2 



R 



k=0 
{k+l)h 

kh 



{m — k — l)h 



sup \{{F'{Xkh)n-A)-'^^w, 

\w\\[f<l 

'm-2 ^(fc+i)/i 



ds 



I fc^o ^fc/x 

+ \\e^(^-'^>^)Okh-OkhhHn;H)fds 



1/2 



and 



,6(l + log(M)) 



m-2 



A;=0 



(m — k — l)h 



{k+l)h 
kh 



sup ||(F'(XfcO)*(-A)-i/2^||^^^_^^^^^^ 

\w\\h<1 



X e 



A(s-fc/i) 



fc/l|lz)((_A)-l/2) 



L2(n;M) 



E 

fc=0 



2 /.(fc+l)/i 



(||Os - Ofch,||L2(n;l/) 



1/2 



+ e 



A{s-kh) 



I)Okh\\L2(n;H)) ds 



for every m G {0, 1, . . . , M} and every M G N. Using now condition (13) in 
Assumption 2 shows 



\\Xmh - Ym \\l'^{Q.;H) 

< AR 



6(l+log(Af)) 



M2^ 
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m-2 



(m — k — l)h 



(31) 



(fc+i)fe 

||c(l + ||^fc/i|lD((-A)i/2)) 



Ikh 

X IKe'^^'"''''^ -^)Ofch|lD((-A)-V2)|lL2(f,.K)ds 

{"^-2 /•{k+l)h 
{R{s-khf 
u—n J kh 



k=0 ' 



1/2 



+ (s - A;/i)'^||Ofc/,||i2(j^.j)((_^)7)))^ds 



and therefore 

\\Xmh - \\L'2{n;H) 

m-2 |-(fc+l)?i 



X IKe^^'"''''^ --^)OM|lL4(j,.^((_^)-i/2))ds 

■ _n ^fc/x I 



fc=0 

for every m € {0, 1, . . . , M} and every M G N. Hence, we obtain 

\\Xmh — \\L'2{n;H) 

,6(l + log(M)) 



<4i?' 



m-2 



k=0 ^ ' 



r{k+l)h 

X / \\ie^^'-"'^-I)OkhhHn;Dii-Ar^n)^' 

J kh 

m-2 .(k+l)h ^ 



fc=0 -^^^ 
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and 



W^mh - Ym \\l2{Q-H) 

< m 



6(l+log(M)) 



cR{l + R) 



(m — k — l)h 

k=o ^ ' 



{k+l)h 



Ikh 

X \\Okh\\Li{Q;D{(-A)^)) ds 

for every m G {0, 1, . . . , M} and every M G N. This yields 

W^mh - ym\\L2{n;H) 



(l + log(M)) 
MS'' 



and hence 

W^mh - Y^\\L^{n;H) 

for every m G {0, 1, . . . , M} and every M G N. Therefore, we have 

™-2 2i?4 r{k+l)h / \ (7-1/2) 
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and, finally, 



W^mh - Y'm\\L2{Q;H) < 6i?' 



6(l + log(M)) 

M2« 



m~2 



k=0 



(m — k — 1) 



< 6R' 



6a+iog(M)) 



(32) 



6(l + log(M)) 
M2« 



M 



■ ds 



for every m G {0, 1, . . . , M} and every M € N. 

6.1.2. Spatial discretization error. Due to (28), we obtain 

e^"^"(e-P7v(0) 



''m— 1 



+ ^ (e^^'"''-^^) - P7ve^(™''-^^))F(Xfc/,) 



+ Omh - PN{Omh) 



< \\e^'^\^-PNm\LHn;H) + llO^h -^'7v(0^,)||i.(f,,^) 



+ 



m— 1 



.k=0 



L2(Q;_ff) 



and 



< lie - PN{0\\L^{n;H) + l|Om/i " (Omh.) ||L2(n;/f) 



'' m— 1 



+ ^Ile^^'"''-'^'^) -P^e^('"'^-'=^)||^(^)||F(X,,)||^.(^^^) 



vA;=0 



< IIU - PNmL^n;H) + ^7v)0,nh||L2(Q;H) 



^ m—1 



Rhi J^||(7-P^)e 



A(mh-kh) i 



lL(H) 
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and hence 

W^m - PN{ym)\\L-^{n;H) 

<\\{-A)-^{I - PM)\\L(^H)\\{-AViU^^n,H) 

+ \\[-A)-^{I - PM)\\LiH)\\{-AVOmh\\LHn,H) 

(m-l \ 
k=0 / 

< {XN)-^{\\{-AraLHn;H) + \\{-Ar0^h\\LHn;H)) 

V k=0 J 

for every m € {0, 1, . . . , M} and every M G N. Therefore, we have 

\\Y^ - PNiym)\\L2{n;H) 

< i^Ny^ i\\^\\L'^(n;D{(-Ap)) + \\Omh\\L'^{n;D({-Ap))) 

(m—l ^ ^ 



and 



\\Y^ - PNiY^)\\LHn;H) 



(33) =i?(A^)-7(^2 + /i(i-^)(l+[^— 

< r{Xn)-^ [2 + (Y^:^ ) - 3^ (^^)"^ 

for every m E {0, 1, . . . , M} and every M G N. 
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6.1.3. Lipschitz estimates. Note that 1^'^^ : —)• F satisfies 



mh ) 



>fc=0 



for every m £ {0, 1,. . . ,M} and every N,M gN. Indeed, in the case m = 
we have 

= e^'iPMiO) + /if j;P^e^(°-^'^)F(y,^'*^)') + P^(Oo) 
\fc=o / 

for every A^, M S N. Moreover, if (34) holds for one m £ {0, 1, . . . , M — 1}, 
then we obtain 



= e^^Y^'^ + h ■ P^e^'^FiY^^''^) + PN{0^m+i)h) " e^'^'^(On^h) 



= e^'^iY^^''^ - PNiO^h)) + h ■ PNe^^F{YZ'''^ ) + P^(O(^+i)0 
and 

+ h . PNe^'^FiY^'^) + PN{0^m+l)h) 

^ P^e^«'"+i)^-'='^)F(y,^'^^) 

k=0 ) 
+ h ■ PNe^'^FiY^'^') + PN{0^m+l)h) 



\k=0 



+ PNiO(^rn+l)h) 



for every N,M gN, which shows (34) by induction. In the next step, (34) 
yields 



/m-1 N 

= hiY,PNe^^""'-'"'^F{X,h) 



, fc=o 
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/ 771—1 

Kk=0 



Vfc=0 / 

for every m G {0, 1, . . . , M} and every A^, M G N. Therefore, we obtain 

II D (V^\ vN,M\\ 
m— 1 



k=0 

777—1 



(35) < hY,{\\PNe'"^"''-'''^h^j,)\\F{XkH) - '*')llL^(f7;H)) 

fc=0 

777—1 



k=0 

777—1 



for every m G {0, 1, . . . , M} and every A^, M G N. Combining (32), (33) and (35) 
finally yields 

II -V- taTV.Mm 

W-^mh — J^rn \\L'^{p.;H) 

< \\Xmh - YmllL'^iniH) 

+ 11^7^ - PN{ym)\\L2{n;H) + l|-PAf(^7^^) " ^7^'*^ II L2 

- M29 7T-y7+c/i2^ ll^fc/7->fc IIl2(Q;H) 

^ k=0 

for every m G {0, 1, . . . , A/} and every A^, M G N. Hence, Gronwall's lemma 
yields 

(l + log(M)) , , . 1 \ 



II^Tnh - -^777 \\L2(n;H) < J^g + J 



e 



,6(l+log(M)) , 1 



(36) <[7R^ ^ +^^(A^J- 

for every m G {0, 1, . . . , M} and every A^, M G N, which shows the assertion. 
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6.2. Properties of the SPDE (1). 
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Proof of Lemma 1. A standard application of Banach's fix point the- 
orem (see, e.g., Section 7.1 in [2]) yields the existence of a unique adapted 
stochastic process X : [0,T] x Q ^ V with continuous sample paths which 
fulfills (14). Moreover, we have 
ft 



(37) 



F{Xs{u;))dsGDi{-Ay) 



for every t G [0,T] and every a; G fi, since 

[\\(-Are^('-^^F{Xs{u))\\^ds 
Jo 



< 



Jo 



7pA(t-s)| 



\L{H) 



WFiXsiu) 



\H 



ds 



< / it-s)-^\\FiXsiuj))\\Hds 



< / it-s)-^c\\Xsiuj)\\H + \\FiO)\\H)ds 
Jo 

< ( [\-''ds] (c( sup \\Xs{uj)\\h) + \\Fm\H 

\Jo J ^ ^0<s<T ' 



and 



< 



Aye^^'-'^F{X,(u,))\\^ds 

s(i-7) 



< 



(1-7) 

j>(l-7) 

0^ 



c sup \\X,{Lo)\\u)^\\Fm\n 
ci sup ||X,(w)||h) + ||F(0)||//) <oo 



holds for every t G [0,r] and every (x) G $7. Assumptions 3, 4 and (37) hence 
imply XiiuS) G D{fy—A)'^) for every t G [0,T] and every a; G Furthermore, 
we have 



-AyXtWu 



< 



-^)^e^*elk+ f\\i-Are^^'~''^FiXs)\\Hds + \\{-ArOt\\H 
Jo 



<\\i-ArC\\H+ r||(-^)V(*-)|L(^)||F(X,)||Hds 
Jo 



10 

+ sup \\i-AyOs\\H 

0<s<T 
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<\\{-Ay^\\H+ [\t-s)-^\\F{Xs)\\Hds+ sup \\{-AyOs\\H 

Jo 0<s<T 

<(\\{-Arc\\H+ sup wi-AyOsWH) 

+ f\t-s)-^{c\\Xs\\H + \\Fm\H)ds 
Jo 

for every t £ [0,T]. This yields 

\\{-Ayx,\\H 

<c [ {t-syWXsWHds 
Jo 

+ (\\i-A)\\\H+ sup Wi-ArOsWH + WFmni fs-^ds)) 

\ 0<s<T \Jo J J 

<(\\i-Arc\\H+ sup ii(-A)^o.ii^+ ^^'7?^?^"^ 

V 0<s<T (I-7) 

+ c\\i-A)-^h^H) [\t-s)-^\\{-Arx,\\Hds 
Jo 

for every t G [0,T]. Hence, Lemma 7.1.1 in [8] shows 

sup wi-AyxtWH 

0<t<T 

< (r(c|| {-A)-^u^H)r{i - 7))'/('^^) ) 

-AraH+ sup ii(-a)^o.iih + ^^7?^?^"^ 



and therefore 

-AVXAIh 

L'1(Q;M) 



sup \\{-AyXt\\H 

0<t<T 



< 



i?(i_,)(r(c||(-A)-^||i(^)r(i-7))'/('-^)) 



X 



-AraLHn;H)+ sup ||(-A)^0,||^^ 

0<s<T 



r(i-^)||F(o)||^^ 



(1-7) 

which shows the assertion. Here : [0,oo) — )• [0,oo) is given by 



< 00, 



^ ^("(1-7)) 



„=orKi-7) + i) 

for every x £ [0, 00) where F : (0, 00) — )• (0, 00) is the Gamma function. □ 
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Lemma 2. Let Assumptions 1~4 be fulfilled. Then we have 

L(H) 



\e — e T.m\ S 



h 

for every ti,t2G (0,T] with ii < ^2- 
Proof. By definition, we have 

= ||(^(t2-ti))-'(e^(*^-*^)-/)IL(^) 

V II Ai P^*i|l it2-tl) 

X ll-4iie IIl(h) 

< sup ^ 1 ( sup xe ' 

\xe(o,oo) ^ / ^ze(o,oo) 

< it2-h) 



h 

for every ti,t2 G (0,r] witli ti < t2- D 

Lemma 3. Let Assumptions 1-4 be fulfilled. Then we obtain 

sup 77 TTe < 

0<ti<t2<T {t2-tl) 

where 6 G (0, ^] is given in Assumption 3 and where X : x [0, T] — t- D{{ 
is the solution of the SPDE (14)- 

Proof. First, let R G [0,oo) be tlie real number given by 

R-=U\\L^in;D(A))+ sup ||F(Xt)||L2(n;//) 

te[o,T] 
( \\0t2 - \\L'2(n-H) 

+ sup I 

0<ti<t2<TV {t2-tl) 

which is finite due to Assumptions 1-4. Then we have 

(38) = ||e^*i(e^(*-*^)e - OWlHu^h) < lle^^*^"*^)? - C\\LHn;H) 
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for every < ti < t2 ^ Moreover, we obtain 

e^(*2-^)F(X,) ds - r e^^''-'^F{X,) ds 
Jo Jo 

F{Xs)ds+ f\e 
Jo 



)FiX,)ds 



< 



and hence 



\\HXs)\\LHn;H) ds + 



)FiXs)ds 



L'2{n;H) 



I ' e^(*^^^)F(X,) ds- r e^(*i-^)F(X,) ds 
Jo Jo 

<R{t2-h)+ f\\e^^'' 
Jo 

<R{t2-ti) + R ['\\e 
Jo 



^\L(^H)\\F{Xs)\\L^in;H)ds 



A{t2s) _ A(ti-s)ua-9)u A{t2-s) _A(ti 



\L(H) 



\L{H) 



ds 



for every < ti < t2 ^ This yields 



e^^^^-'^F{Xs) ds - r e^^^'-'^F{Xs) ds 
Jo Jo 

< R{t2 -ti) + 2R{t2 -h)^ ih - s)-^ ds 

Jo 



due to Lemma 2 and therefore, we obtain 



t2 



oMt2-s) 



FiXs) ds 



ti 



(39) 



FiX,)ds 



L'2{n;H) 



ft! 

< R{t2 -ti) + 2R{t2 -hf / s"^ ds 

Jo 

< R{t2 - h) + (^^iir(i-^)) {t2 - hf 



for every <ti < t2 <T. Combining (38), (39) and Assumption 3 yields the 
assertion. □ 

Lemma 4. Let Assumptions 1~4 be fulfilled. Then we obtain 
\\{Xt2 -Ot^)-{Xt^ - OtMb^in^H) 



sup 

0<ti<t2<T 



{t2-tl) 



< oo, 
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where O:[0,T] x — ^ D{{—A)"') is given in Assumption 3 and where X : 
[0,T] D{{-Ay) is the solution of the SPDE (I4). 

Proof. First, let R G [0,oo) be the real number given by 

R-=U\\L^n;D{A))+ \\^(.^t)\\mn;H)+ SUp , ' , 

ie[0,T] 0<ti<t2<T (i2-cij 

which exists due to Lemma 3. Then we have 

(40) <\\e^'-''-''k-aLHn;H) 

<R{t2-ti) 
for every <ti <t2 <T. Moreover, we have 



f ' e'^^^^-'^FiXs) ds- r e^(*i-^)F(X,) ds 
Jo Jo 

e^(t^-')F{Xs)ds 

Jti 



L'2{n;H) 



+ 1 (e^(*2-^) 



L2{n-H) 



rt2 


L 


||g^(t2-s) 










+ 






Jo 






+ 






Jo 



■,H) 



ds 



Ait2~s) _ e^(<i-^))(i7(x,) - F{Xt,))ds 



and therefore 



f ' e^(*2-^)F(X,) ds - e^^''-'^F{Xs) ds 
Jo Jo 

<R{t2-ti) 



\L'2(n;H) 



ds 



+ R 



"(e^(*2-«)_e^(*i-^))ds 



LiH) 
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for every <ti <t2 <T. Hence, we obtain 



t2 



F{Xs)ds 



FiXs)ds 



L2(n;_ff) 



<R{t2-ti) 

+ C 







ds 



and 



R 



FiX,)ds 



FiX,)ds 



L{H) 



<R{t2-ti) + cR [ \\e 
Jo 



A(t2-s) _ pA(ti-s) 



\LiH) 



\s-tifds 



+ R 



t2 



e^'ds 



'{t2~h) 

for every < ti < t2 ^ This shows 



t2 



,A(t2-s) 



FiXs)ds 



,A{ti-s) 



e^'ds 



FiXs)ds 



L(H) 



L^{n;H) 



< R{t2 - h) 



s-tif ds + R 



t2 



e^' ds 



ih - s) 
nd 

r e^(*2-^)F(X,) ds - r e^(*i-")F(X,) ds 
Jo Jo 



(t2"tl) 



e^' ds 



L{H) 



due to Lemma 2 and 

rt2 



fti 

< R{t2 - ti) + cR{t2 - ti) / {ti - ds + 2R{t2 - ti) 

Jo 

r-ti 

= R{t2-h) + cR{t2-ti) / s^'^~^Us + 2R{t2-h) 

Jo 

<{R + cR{T + 1)9^^ + 2R){t2 - h) 
for every < ti < t2 ^ Combining this and (40) shows the assertion. □ 
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